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We present a theory of phase transition in quantum critical paraelectrics in the presence of a 
quenched random-Tc disorder using replica trick. The effects of locally ordered regions and their 
slow dynamics are included by breaking the replica symmetry at a vector level. The occurrence of a 
mixed phase at any finite value of disorder strength is argued and a a broad power law distribution 
of quantum critical points is predicted. Results are interesting in the context of a certain class of 
disordered materials near quantum phase transitions. 
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I. INTRODUCTION 



Following the recent progress in the understanding 
of quantum phase transitions in ferroelectrics^ and 



the strongly correlated systems with disorder——, we 
make an attempt to explain certain features of dis- 
ordered quantum paraelectrics. In general, a dilute 
amount of quenched impurity and associated disor- 
der, can create locally ordered regions even above 
the transition point of the corresponding pure sys- 
tem. Near a phase transition, large size droplets 
become more probable and their slow dynamics be- 
comes an important factor to determine the nature 
of quantum phase transition in the disordered sys- 
tem. Some dielectric materials, like SrTiOa, KTaOa 
do not show any ferroelectric transition even at 
zero temperature in their pure form. These rather 
end up with a very high, temperature independent, 
static dielectric susceptibility at low temperatures^. 
The absence of a transition and the saturation at a 
high dielectric susceptibility in these materials cor- 
respond to the smallness of the zero temperature 
gap in an optical branch and are believed to be near 
a ferroelectric quantum critical poin^. These ma- 
terials have been studied in the past several years 
and the role of quantum fluctuations in the sup- 
pression of ferroelectric transitions and the efi^ects 
of pressure^ and impurity^ were emphasized in ear- 
lier literature. A relevant question in this line would 
be the effects of impurity and disorder in these sys- 
tems. Most of the disordered quantum paraelectrics 
show relaxor behavior which is often described as a 
classical glassy behavior of a dipolar system. In this 
work we focus on the effects of disorder in quantum 
critical behavior of certain ferroelectrics. Such issues 
were addressed in case of classical critical behavior 
earlier— and some attempts to make a quantum gen- 
eralization of it in the context of itinerant magnets 
have been proposed in the recent pasb^"— . We use 
some of the earlier results and develop a new mean 



field description of the possible low temperature be- 
havior of a disordered quantum critical paraelectric. 



II. THEORY 

Polarization fluctuations in these materials are well 
described by a zone center transverse optic mode. 
In a Heisenbcrg spin(dipole) model with dipolar in- 
teractions, the longitudinal modes are gapped out 
and transverse fluctuations arc able to become soft. 
For simplicity, we restrict ourselves to one transverse 
branch. Assuming that the system is soft in one 
transverse direction, we start with a one component 
Landau-Ginzburg- Wilson quadratic action, suitable 
for near quantum critical paraelectrics. 

Apure = -^^(tJ?, + q^ + r)|(/i(a;„,q)p 
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Here /3 = l/ksT where fc^ is the Boltzmann con- 
stant and T is the temperature. Here a;„ = 2'KnT 
with integer n, is the Bosonic Matsubara frequency. 
The parameter r determines the gap in polarization 
fluctuations and r = is the mean field quantum 
critical point of the pure system. Disorder is intro- 
duced into the problem as a random variation of r 
in real space and the disorder contribution to the 
above mentioned quantum action is given by. 
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where (5r(x) is assumed to follow a Gaussian proba- 
bility distribution with variance viz.. 



P((5r(x)) cx exp{—— j (ix(5r^(x)}. 



(3) 
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so that (Sr(x) = and Sr{x)Sr{y) = g6'^{x — y). To 
calculate the disorder averaged free energy, we intro- 
duce replicas of the order parameter (pa with replica 
index a = l,....,n and write F = — — l)/n, 
taking n — > at the end of the calculation. Dis- 
order averaging generates interactions between fields 
with different replica indices and the resulting action 
becomes, 
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Here a, b arc the replica indices which take positive 
integer values up-to some integer n. 
To begin with, first we consider a replica symmetric 
case. We define replica symmetric solution as replica 
independent field configurations, i.e. 



(f>a (x, r) = '?5'(x, r) for all a 



(5) 



and their replica diagonal two point correlation func- 
tions, i.e. 



Xab = XaaSab = Xof^ab for all a, 6. 



(6) 



We consider a paraelectric phase i.e. < >= 
and make a self-consistent quasi-harmonic approx- 
imation to decouple the quartic term. In this 
scheme a quartic term such as J dxdT(j)'^{x,T) is 
decomposed as Ao / dxdrcjp {x,t)). Where Aq = 
/ dx dT'4>^{x ,t'). Thus susceptibility of the disor- 
dered paraelectric can be written as. 
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(w^ + q2 + r + Ao) ■ 



(7) 



In the above equation Ao describes the fiuctuation 
corrections to ferroelectric gap and is defined by the 
following self-consistent equation 

Ao = ^{uxo{^m,q) - 9XoiO,q)) 

m,q 

= uj d3gi_coth/3^],-g J d^g-L, (8) 

where the fluctuation corrected natural frequency flq 
is dcflned as, 



nl = q^ + r + \o. 



(9) 



It is to be noted that the second term in equation 
[5] is a zero frequency contribution. The reason is 



that we consider quenched disorder which has no dy- 
namics and thus strongly correlated in time. How- 
ever the above two equations can be obtained by 
integrating 6r{x) without introducing replica trick 
and need to be solved sclf-consistently. It is clear 
form the expression for Aq that the second integral 
in the equation ([8|) gives a shift in gap and de- 
pending on its strength controls quantum fluctua- 
tions. In this scheme the solution of the equation 
r — g J d'^qfl~^ = for r gives the quantum critical 
point. Right at that point various physical quanti- 
ties follow power law dependencies in temperature, 
e.g. static dielectric susceptibflity ~ T"^. The ef- 
fects of disorder considered here are identical to the 
effects of hydrostatic pressure as discussed in our 
previous worki. 

To include the effects of spatial inhomogcneity cre- 
ated by disorder we need to break the replica sym- 
metry at the vector leveUS. In this approximation 
field configurations are assumed as 

(?!)a(x, r) = 0fe(x, t) + ?/'(a;, T)fora = 1, 
(f)a{x,T) — iP{x,t) a ^ k + 1, , n, (10) 

and the correlation functions to be block-diagonal 

Xabix,T) ^ xi{x,'r)+X2ix,T)SabioTa,b^l,..,k 
= X2{x,T)Sabioia^ k + 1, ,n. (11) 

Here /c > 1 is an integer that determines the de- 
gree of the symmetry breaking process. In classical 
treatment for gk > u, (f>i; corresponds to localized 
solution given as. 



-d'{^/rx). 



(12) 



p{gk - u) 

So that Tpiz) obeys a scale independent equation 

-V^V(^) + ^(^)-^'(^)-0. (13) 

The proper boundary conditions are V'(O) ~ con- 
stant and tJj{±oo) = 0. Equation 1131 has exponential 
decaying solutions for x 3> \/r and is smooth for 
x < y/r. The size of the droplet R is determined 
by the dipolar correlation length and R ^ At 
very low temperature the dynamics of the droplets 
become important. In a simplest approximation spa- 
tial and the time dependent parts of the polarization 
field can be decoupled completelj^. 



{x,t) = <pk{x)T(T). 



(14) 



Within this approximation the dynamics of the lo- 
calized solution can be cast in a problem of an un- 
damped Bosonic particle in a double well potential 
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with frequency = 2 in some dimensionless unit. 
Authors of referencei^ derived a tunnel splitting of 
ground state energy 



(15) 



with ro ^ E1/E2 a constant, where En = 
J dz(j){z)'^^ . Integrating out the fluctuations due to 
droplets in a Gaussian approximation, an effective 
action for the paraelectric fluctuations "0(0;, r) can 
be written as, 

m,q,a,6 

The presence of droplets introduces a "gap-matrix" 
{M.ab] which contains k y. k block with elements, 

9k — 3u^.^ 2qkr 

Mab--r{l-'^ \L)Sab-^ Ai (17) 

gk — u gk — u 

and diagonal elements for the remaining n~ k rcpli- 



gk 



gk — u 



^L)5ab 
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Here \l encodes the contributions from the localized 
solutions along with their dynamics and is given as 

\l = dz< iP{z)T{io)^{z)r{io) > 



= I dzr{z)^<T{io)T{u:)> 



— at T = 



(19) 



Here a;± = 2 ± r^. It is to be noted that the 
vector breaking of replica symmetry not only in- 
troduces inhomogeneous solutions but also glassy 
effects through off-diagonal elements in the gap- 
matrix. Putting Al = identically, we get back the 
behavior of a pure system. However in this scheme 
replica correlators for disordered paraelectric is given 

by 

Xlbi^^n, q) = + Ci^)5ab + Mab)- (20) 

We look for replica diagonal correlations in equations 
pT)) and (flSl) . Diagonalization of the gap- matrix are 
given as, 



(l-^A,),a = l,..fc-l, 
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Using equation ([5T]) and ([T^ we find the values of 
r at which zero temperature diagonal susceptibility 
(~ 'j^) diverges. The instability point depends on 
the disorder strength and the value of k and is given 

as, 



-ro/log(l 
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(22) 



fc + 1 , . . , n 



and is k dependent. Here A is a, system depen- 
dent parameter. For a simple minded analysis, let 
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FIG. 1: Schematic phase diagram at zero temperature 
in g — r plane for k = 1,2 and 4. Above a certain value 
of g there is an ordered phase for r > . Moreover the 
region between lines corresponding to various values of 
k represents a mixed phase. 

us consider the a ~ k + I, ..,n elements only. For 
I (g^^^^) I << 1, can be written in the following 
form 



rc(fc) = To 



{gk - u) 
gk 



= B-^{B = r,,C^'-^). (23) 
K g 

Since the choice of k is random, depending on its 
distribution at the limit n — > 0, we can estimate 
a distribution, hence width of Tc- In a replicated 
action with n replicas, k can be chosen in ways. 
Thus wc can define a normalized distribution of'P{k) 
as follows 



V{k) 



cv 



1 



r(n) 



(24) 



ELi^r 2"-ir(fc)r(n-A:)- 

Since the gamma function with negative argument 
is infinity, the limit of k can be extended to infinity. 
In the limit n — ^ 0, using the asymptotic form of 
gamma functions, 'P(fc) can be approximated as^^^ 



V{k) 



1 (-1) 
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1 cos TT COS Trfc 
log 2 k ' 



(25) 
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Negative values of 'P{k) for some values of k may 
turn out to be counter intuitive to the usual notion 
of a distribution function. But such distributions are 
allowed in replica scheme. There are several possible 
broken replica symmetric cases, each characterized 
by the number k which follows a distribution V{k). 
For a fixed disorder strength 5, each k results a dif- 
ferent instability point re- Instead of k, if we char- 
acterize various possible broken replica symmetric 
cases by Tc, a distribution of Tc can be estimated as 



/\k 1 



X cos(7rA:). (26) 



This is a broad power-law distribution of around 
a system dependent parameter B with a cosine fac- 
tor. The probability distribution can be assumed to 
be smooth around k = any positive integer, exclud- 
ing zero. The expansion around fc = is excluded 
as it corresponds to small u/g limit where the ac- 
tion |4] becomes unstable even in a replica symmetric 
ansatz. In that limit the system will undergo a first 
order transition in a replica symmetric analysis, the 
stability of the system needs a 0^ term in the action 
Ulwhich will lead to more complicated localized solu- 
tions in a broken replica symmetry picture. However 
we focus on those u/g values where the above possi- 
bilities are not present and the distribution function 
is smooth. It is to be noted that the power law na- 
ture of ■P(rc) arises because of the dynamics of the 
locally ordered regimes and also depends on the dis- 
tribution of k used. Neglecting cosine factor within 
some range of Tc say {B + R, B — R) , average suscep- 
tibility of the disordered quantum paraelectric can 
be estimated as, 



x{r,T) 



B+R 



dr. 
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(27) 



It is evident that inclusion of fluctuations due to 
locally ordered regime introduce a parameter R ^ 
0{u/g) and changes the usual quantum critical be- 
havior of a paraelectric. In the limit r ^ B the 
temperature dependence of a disordered quantum 
paraelectric can be predicted as, 



— — — ■ J constant, T « R 
'^\l/T\ T»R. 



(28) 



This is a deviation from the standard quantum crit- 
ical behavior which predicts x{T) ^ in a mean 
field theory. 

III. DISCUSSIONS 



In this work, the low temperature dielectric behavior 
of a quantum paraelectric in presence of quenched 
disorder is addressed. A suitable action for these 
materials, with random Tc type disorder have been 
studied using a replica trick. The effects of disor- 
der induced locally ordered regimes and their tun- 
nelling in the low temperature are captured in this 
formalism. We derive an expression for the distribu- 
tion of instability points for a fixed value of disorder 
strength and demonstrate the possibility of a mixed 
phase at non-zero disorder strength. This analy- 
sis shows that the instability points follows a broad 
power law distribution around a system dependent 
parameter with a cosine correction. By using such 
distribution we are able to show analytically how the 
temperature dependence of static dielectric suscep- 
tibility of a disordered quantum critical paraelectric 
deviates from its pure counterpart. Our analysis is a 
completely new attempt in the context of the effects 
of disorder in ferroelectrics near a quantum critical 
point. In a qualitative manner it predicts certain 
new features such as occurrence of a phase with mix- 
ture of critical and non-critical regimes with a dis- 
tributions of transition points which are missing in 
earlier works in similar issues in context of itinerant 
magnets. Moreover the whole analysis is interesting 
in context of the use of replica trick to incorporate 
disorder induced inhomogeneities or locally ordered 
regime in the studies of quantum phase transition 
and may turn out to be useful in explaining cer- 
tain experimental results on disordered ferroelectrics 
near a quantum critical point. 
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